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ABSTRACT 

We examine 11 XMM-Newton observations of the giant spiral galaxy NGC 1961, 
allowing us to study the hot gaseous halo of a spiral galaxy in unprecedented detail. 
We perform a spatial and a spectral analysis; with the former, the hot halo is detected 
to at least 80 kpc and with the latter its properties can be measured in detail up to 42 
kpc. We find evidence for a negative gradient in the temperature profile as is common 
for elliptical galaxies. We measure a rough metallicity profile, which is consistent with 
being flat aX Z ^ 0.2 Zq. Converting to this metallicity, the deprojected density prohle 
is consistent with previous parametric fits, with no evidence for a break within 42 kpc 
(^O.lRvir)- Extrapolating to the virial radius, we infer a hot halo mass comparable 
to the stellar mass of the galaxy, and a baryon fraction from the stars and hot gas of 
around 30%. The cooling time of the hot gas is orders of magnitude longer than the 
dynamical time, making the hot halo stable against cooling instabilities, and we argue 
that an extended stream of neutral Hydrogen seen to the NW of this galaxy is instead 
likely due to accretion from the intergalactic medium. The low metallicity of the hot 
halo suggests it too was likely accreted. We compare the hot halo of NGC 1961 to hot 
halos around isolated elliptical galaxies, and show that the total mass determines the 
hot halo properties better than the stellar mass. 

Key words: galaxies: haloes, galaxies: spiral, galaxies: individual: NGC 1961, X-rays: 
galaxies 


1 INTRODUCTION 

X-ray emission from hot gas appears to be a generic fea¬ 
ture of massive halos. Hot gaseous halos suffuse nearly every 


galaxy cluster and many galaxy groups (Forman & Jones 


1982 


|Sarazin||1986[ [Sun et al.||2009| [Kravtsov fc Borgani 


20121, and they are also very common (possibly ubiquitous) 
around massive elliptical galaxies I Forman et al. 


biano|19'89 1, including field ellipticals (Anderson et al. 20151. 
The hot halo is strongly affected by both feedback from the 
galaxy and by accretion from the intergalactic medium (e.g. 


Cen & Ostriker 2006 Roncarelli et al. 20121. These pro¬ 


cesses, which are fundamental for understanding galaxy for¬ 
mation, can therefore be studied through X-ray observations 
of the hot halos. 

In the aggregate, X-ray emission can be described by 
power-law relations as a function of stellar mass (e.g. [Hels- 


Reiprich & Bohringer 20021, although the scatter in these 


relations is considerable. Potentially active galactic nucleus 
(AGN) feedback (e.g. Churazov et al.|2000 Churazov et al. 
20011 might also be important. For example, the slopes of 


these power-law relations differ from the self-similar predic¬ 
tion, which can be used as a constraint on the effect of AGN 


feedback (Gaspari et al. 2014 Anderson et al. 20151. For 


clusters and groups the luminosity and temperature of the 
gas are also related (e.g. [Mitchell et al.[1977| Mushotzky et 
al.||1978[ [David et al.jfTAsi [Bryan fc Norman||1998| [Pratt 


et al. 20091, but this relation seems to break down at the 


scale of galactic halos (jFukazawa et al.||2006[ [Humphrey et 


al.|2006| Diehl fc Statler|2008 1. Some of the scatter in these 

relations may be correlated with the degree of rotational 
support in the central galaxy as well ( [Sarzi et al.|2013l [Kim[ 
fc Fabbiano[[2013 (. 


don et al. 

2001| [O’Sullivan et al. [[20031 [Mulchaey & Jel- 

Another important issue is the baryon fraction of mas- 

tema|2010 

Boroson et al.|2011 

1 and halo mass ( 

Kaiser] 1986 

sive halos. In galaxy clusters, the hot halos contain the 


overwhelming majority of the baryons associated with the 


email: michevan@mpa-garching.mpg.de 


cluster (Ettori 

2003 Gonzalez et al.||2007 Dai et al. 2010 

Lagana et al. 

20131, although the stellar phase begins to 
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2 Anderson, Churazov, and Bregman 


grow in importance as the halo mass decreases. An open 
question is the contribution of the gas phase in group and 
galaxy halos. If it is sufficiently massive, then the baryon 
fraction remains constant from galaxy clusters (which are 
nearly baryon-complete; [White et al.|19^ l down to isolated 
galaxies. Current censuses appear to show a declining baryon 
(stars + ISM + hot gas) fraction as the halo mass decreases 


ijMcGaugh et al. 2010| Anderson fc Bregman 2010 Papaster- 


gis et al.||2012 1, suggesting a problem of “missing baryons’ 

from galaxies, but significant questions remain. 


Studies by Planck Collaboration et al. (20131 and Greco 


et al. (20141 find a self-similar relation as a function of halo 


mass for the hot gas pressure as probed by the Sunyaev- 
Zel’dovich (SZ) effect. This suggests that, at least within the 
Planck beam (typically several galactic virial radii), the Gos- 
mic baryon fraction is approximately recovered in hot gas 
(although see Le Brun et al.|2015 l. There are also claims of 
individual elliptical galaxies with extremely massive gaseous 
halos that bring the baryon fractions of the systems up to 
the cosmological value ( Humphrey et al.|2011 Humphrey et 
al. 20121, although for at least one of these claims it has been 


shown that the spectral modeling employed for the analysis 
violates X-ray surface brightness constraints and therefore 


signihcantly overestimates the total hot gas mass (Anderson 
fc Bregman||2014 1. Many simulations predict that the den¬ 
sity profile of the hot halo should systematically flatten as 
the halo mass decreases, due to the increasing importance of 
galactic feedback which inflates the entropy of the hot halo 
( White fc Frenk[|1991 Benson et al.|2000 l. A flatter density 
profile would have less pronounced X-ray emission, poten¬ 
tially reconciling the SZ results with the X-ray constraints. 

Around spiral galaxies, the picture is a bit different. 
Starbursting galaxies produce X-ray emitting winds ([Strick¬ 


land fc Stevens|2000[ [Strickland et al. 12004] |Tullmann et~ 


2006 Li & Wang 20131, but extended gaseous halos have 


generally proven to be very difficult to detect around spiral 


galaxies (Bregman & Glassgold 1982 Benson et al. 2000 


Rasmussen et al.|2009 Bogdan et al.|20T5 I. The exceptions 

are the most massive spiral galaxies. Extensive hot halos 
have been detected around the giant spirals NGG 1961, UGG 


12591, NGG 6753, and 2MASX J23453268-0449256 (Ander- 


son^_&_^regmmij[20^ [Dai et aL]|2012| [Bogdan et ai.||2013[ 


Walker et al. 20151. There is also a weak ROSAT detec¬ 
tion of hot gas around NGG 266 ( Bogdan et al.|2013b l, and 
suggestion of extended hot gas around a stack of ROSAT 
images of nearby isolated spirals ( Anderson et al.|2013 |. 

It is unclear what explains the difficulty of detecting hot 
gaseous halos around massive spiral galaxies. Weak lensing 
studies (e.g. Sheldon et al.|20'04 Hoekstra et al.|2005 Man- 


delbaum et al.[2006 Velander et al.[2014 1 generally find that 
massive blue centrals and massive red centrals obey different 
M * —Mhaio relations, so that at fixed stellar mass, massive 
blue galaxies lie in less massive halos than red galaxies. If 
the gaseous halo is responding primarily to the potential 
of the dark matter halo, this could plausibly explain the 
relative X-ray faintness of hot halos around spiral galax¬ 
ies. Alternatively, if the hot halo is tightly coupled to the 
galaxy (through feedback, accretion, or some combination 
of the two), then the differences in feedback and accretion 
between spirals and ellipticals might explain the different 
hot halo properties. 

Around the massive spirals, the hot gas generally ap¬ 


pears to be smoothly distributed with approximate az¬ 
imuthal symmetry, and has been detected to radii of 60 kpc 
or so from the galaxy (although there is no indication that 
the halo truncates here). The temperatures of the hot ha¬ 
los are generally close to the expected virial temperatures 
of these systems, suggesting the hot gas is roughly in hy¬ 
drostatic equilibrium with the halo potential. Gomparisons 
with cosmological simulations have been able to reproduce 
the X-ray emission around spirals (Li et al. 2014[ Bogdan 


et al. 20151, although these comparisons have been unsuc¬ 


cessful in reproducing the X-ray properties of spirals and 
ellipticals simultaneously. A unified picture of the X-ray cir- 
cumgalactic medium is still lacking. It is particularly im¬ 
portant to push the X-ray observations of spiral galaxies to 
larger radii, to probe the regime where most of the mass in 
the halo is predicted to lie, and to be able to constrain the 
properties of the dark halo and the hot halo simultaneously. 

The real quantities of interest are the radial density and 
temperature profiles of the gas (or equivalently the entropy 
and pressure prohles), since these encode the feedback his¬ 
tory of the galaxy and determine the behavior of the hot 
halo. The temperature profile can easily be measured from 
the X-ray spectrum, assuming enough photons are available 
to divide the field into several annuli. So far no observations 
have been deep enough to permit this measurement for the 
hot halo of a spiral galaxy. 

The density profile can be inferred from the X-ray sur¬ 
face brightness profile, using an assumed gas temperature 
and metallicity profile and a few assumptions to deproject 
the observations. Typically spherical symmetry and a fiat 
abundance profile are assumed, as well as either a flat tem¬ 
perature profile (leading to a “beta profile” for the density; 
[Cavaliere fc Fusco-Femiano||1978| ) or a temperature profile 
appropriate for a cool-core galaxy cluster (leading to a modi¬ 
fied “beta profile” for the density; Vikhlinin et al.[2006 f The 
assumed metallicity has a significant effect on the final re¬ 
sult, since the plausible range of gas metallicities spans more 
than an order of magnitude (roughly 1/10 Solar to slightly 
super-Solar). For metallicities in this range, the conversion 
from soft-band surface brightness to density depends ap¬ 
proximately on the square root of the gas metallicity. In 
terms of the total gas mass, the extrapolation to large radii 
is even more signihcant since the hot halo is typically only 
detected out to about a tenth of the virial radius; most of 
the inferred gas mass lies at larger radii where we have few 
constraints on the density profile. The farther out the sur¬ 
face brightness profile can be measured, the more reliable 
the extrapolation becomes, leading to better estimates of 
the total hot halo mass. 

In order to improve the constraints on the density pro¬ 
file, and to make a first measurement of a temperature pro¬ 
file, we re-observed the giant spiral galaxy NGG 1961 with 
XMM-Newton for an additional 215 ks (adding to the 74 
ks of observations of this galaxy already taken with XMM- 


Newton and presented in Bogdan et al. (20131). In this paper 


we report our analysis of these observations of NGG 1961. 
In section 2 we discuss the data reduction. In Section 3 we 
present a spatial analysis of this galaxy, and in Section 4 we 
present a spectral analysis. In section 5 we combine these 
analyses and measure the pressure and entropy prohles of 
the hot halo. In section 6 we discuss our results in the con- 
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The Hot Halo Around NGC1961 3 


text of the “missing baryons” problem and in comparison to 
isolated elliptical galaxies. 

NGC 1961 is an extremely massive late-type spiral 
galaxy. Based on its recessional velocity of 3934 km/s listed 
in the NASA Extragalactic Database, and an assumed 
Planck cosmology (Planck Collaboration et al. 20151 with 
Hq = 67.8 km/s/Mpc, we estimate the distance of NGC 
1961 to be 58.0 Mpc, so that one arcminute corresponds 
to 16.6 kpc. The K-band luminosity of this galaxy is then 
5.6xlO^iL0, corresponding to a stellar mass of 3 x IG^Mq 
(assuming a M/L ratio of 0.6, which is the rough expecta¬ 
tion based on its K-band luminosity; Bell fc de Jong|2001 1. 
For the B-V= 0.6 color of this galaxy, a Chabrier (20031 
initial mass function (IMF) gives the same M/L ratio, but 
a Salpeter (19551 mass function gives a M/L ratio of 1.2, 


yielding an even larger stellar mass for this galaxy. How¬ 
ever, this latter M/L ratio is disfavored by [McGaugh fc 


Schombert (20151, who use both population synthesis and 


IMF-independent constraints to determine a nearly univer¬ 
sal M/L ratio of 0.57 in the K-band, close to our value of 0.6. 
The inclination-corrected HI rotation velocity is 340 km/s 
(Haan et al. 2008) at a projected distance of 43 kpc (al¬ 
though it reaches higher values - up to about 450 km/s - 
at smaller radii). The virial radius of this galaxy is approx¬ 
imately 490 kpc, as we discuss in section 5.1. We use the 
coordinates from the NASA Extragalactic Database for the 
center of the galaxy. 


Table 1. XMM-Newton Observations of NGC 1961 


obsid 

Texp 

(ks) 

Tmosi 

(ks) 

7mOS2 

(ks) 

TpN 

(ks) 

0673170101 

37.9 

20.4 

20.5 

15.2 

0673170301 

35.9 

13.1 

16.8 

7.0 

0723180101 

23.7 

18.8 

19.2 

14.0 

0723180201 

22.9 

5.4 

5.9 

2.4 

0723180301 

26.5 

9.5 

10.5 

2.4 

0723180401 

23.7 

3.2 

6.8 

1.6 

0723180501 

22.9 

2.7* 

4.0* 

1.6* 

0723180601 

26.9 

6.6 

8.1 

4.8 

0723180701 

22.0 

7.4 

7.1 

5.3 

0723180801 

22.0 

11.9 

12.4 

8.3 

0723180901 

24.5 

10.0 

9.5 

5.6 

total 

289.0 

106.5 

116.8 

65.6 


List of XMM-Newton observations of NGC 1961. The first column 
is the obsid and the second column is the total duration of the 
observation as listed in the XMM-Newton Science Archive. The 
next three columns show the length of the good time intervals 
(GTIs) for each instrument after using the XMM-ESAS software 
to filter each observation, as described in the text. These obser¬ 
vations were heavily contaminated by flaring; much less than half 
of the total exposure time was useable for analysis. Observation 
0723180501 was so heavily contaminated that it was still not use- 
able even after GTI filtering, so we discard this observation from 
subsequent analysis (and do not include it in the GTI filtered 
total exposure time at the bottom of this Table). 


2 OBSERVATIONS AND DATA REDUCTION 


We obtained nine observations of NGC 1961 with XMM- 
Newton, ranging from 22-27 ks in length each. The aim- 
points of each observation were varied such that they formed 
a 3 X 3 grid, with a typical separation between aimpoints 
of about an arcminute. This layout makes the data reduc¬ 
tion and spectral fitting more cumbersome, but it offers the 
advantages of flattening out the exposure map around the 
galaxy and reducing the effects of vignetting. It also allows 
us to identify and separate various non-astrophysical back¬ 
ground components (soft protons, instrumental background, 
solar wind charge exchange) from the astrophysical back¬ 
grounds like the Galactic halo and the Cosmic X-ray back¬ 
ground. We also re-analyzed the two previous XMM-Newton 
observations of this galaxy which were originally discussed 


Bogdan et al. (20131 (hereafter B13). The list of all 11 


observations is presented in Table 1. 

We reduced the data following the procedure out¬ 
lined in the XMM-Extended Source Analysis Software 


(ESAS) Cookbook (Snowden & Kuntz 20141. We used 


HEASOFT V. 6.16 and XMM-SAS v. 13.5.0 for the data 
reduction and analysis. We first ran the initial process¬ 
ing commands (epchain, epchain withoutoftime=true, 
pn-f liter, emchain, mos-f liter) to produce filtered events 
files for each observation. The data were heavily affected by 
flaring, so this processing was especially important. We ex¬ 
amined each lightcurve manually to verify that the pipeline 
processing was working correctly. We also examined the fil¬ 
tered events files for CCDs in anomalous states and excluded 
them. During this process we found that obsid 0723180501 
was so heavily affected by flaring that it was effectively un¬ 
usable; we therefore discarded this observation for the sub¬ 


sequent analysis. Table 1 shows the effects of the initial pro¬ 
cessing on each event file. 

Next, we identified point sources and produced point 
source masks for each observation. Since the fields of view 
overlap in each observation, we produced broad-band images 
and exposure maps from each events file and merged them 
into a single broad-band image. We then ran the Chandra 
wavdetect algorithm on the merged image with default pa¬ 
rameters in order to identify point sources; at a projected 
radius of 3’ (50 kpc) from the center of the galaxy we esti¬ 
mate our limiting point source flux to be around 7 x 10“^^ 
erg s“^ cm“^ in the soft (0.4-1.25 keV) band, correspond¬ 
ing to a luminosity of 3 x 10®^ erg s“^. This is sufficient 
to resolve the brightest X-ray binaries, but there will be an 
important unresolved component remaining. We verified the 
point source list with manual inspection, then passed the list 
of sources to the region and make-mask routines in order to 
construct masks with the appropriate radii for each obser¬ 
vation (setting the radius to encircle 85% of the expected 
energy for each point source). 

Finally, we used the mos-spectra, mos_back, 
pn-spectra, and pn_back routines to generate broad¬ 
band spectra and to generate images in the 0.4-1.25 keV 
(“soft”) and 2.5-7.0 keV (“hard”) bands. These routines 
also generate background files (spectra and images) with 
the estimated particle background in each observation. We 
bin these spectra such that there are at least 20 counts in 
each energy bin. 

To fit the remaining backgrounds, we extract a spec¬ 
trum from each observation after masking out point sources 
and masking out a circle centered on the galaxy with radius 
8’ in order to exclude any possible emission from the galaxy 
or its hot halo. We bin each of these background spectra as 
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4 Anderson, Churazov, and Bregman 


above, and perform a joint analysis of all 30 spectra using 


XSPEC V. 12.8.2 (Arnaud 19961 and the associated Python 


wrapper PyXSPEC. We analyze the MOS spectra in the 0.3- 
11.0 keV energy band and the PN spectra in the 0.4-6.5 keV 
energy band (at higher energies the PN spectra are affected 
by a number of instrumental lines). We use the abundance 
tables of Anders fc Grevesse| (19891. 

The spectral model is complex and is based on the pro¬ 
posed model outlined in the ESAS cookbook. We use two 
Solar abundance APEC models ( [Smith et al.|[^0ip with 
temperatures of 0.1 and 0.25 keV to model the Local Bub¬ 
ble and the Galactic halo respectively. We model the Cosmic 
X-ray background (CXB) with an absorbed P = 1.44 power- 
law, and we set the absorption to the Galactic value at the 
location of NGC 1961 (8.2 x 10^° cm“^, Dickey & Lockman 
1990 Kalberla et al.|2005 1. These components each have free 


normalization, but we tie the normalizations in each obser¬ 
vation together (with appropriate correction factors for the 
differences in the area of each region, and with a free param¬ 
eter which is allowed to range from 0.9-1.1 for each observa¬ 
tion to allow for uncertainties in the calibration, as suggested 
by Snowden & Kuntz (2014l). We also include a number of 
components in our model which account for instrumental 
or time-variable backgrounds, and are therefore allowed to 
very between observations. First is a zero-width A1 Ka in¬ 
strumental line at 1.49 keV with free normalization for each 
detector and each observation, and a zero-width Si Kq in¬ 
strumental line at 1.75 keV with free normalization for the 
MOS detectors in each observation. Second is a power-law 
for the soft proton background, with free normalization and 
slope for each detector and each observation, although the 
slope is constrained to have an index between 0.1 and 1.4. 
This component uses a diagonal response matrix provided 
with the XMM-SAS software package. The ESAS cookbook 
suggests that one may assume the soft proton background 
in each observation has the same slope for both MOS detec¬ 
tors, but given the potential importance of the soft proton 
background due to the significant flaring in our observa¬ 
tions, we choose to relax this assumption and fit the proton 
background for each instrument separately. Finally, we also 
include six zero-width line features to account for the So¬ 
lar wind charge exchange (SWCX) background. These lines 
have energies fixed at 0.46, 0.57, 0.65, 0.81, 0.92, and 1.35 
keV corresponding to Cvi, Ovii, Oviii, Oviii, Neix, and 
Mg XI transitions respectively. The SWCX background is al¬ 
lowed to vary between observations but its normalization is 
fixed across all three instruments during each observation. 

In order to reduce the dependence of the results on ini¬ 
tial guesses, we used the steppar command for each of the 
SWCX lines to explore different values for the normalization 
and improve the fits. The final fit had a reduced of 1.0134 
for 8615 degrees of freedom, which is an acceptable fit. We 
use the normalizations of the CXB, Local Bubble, Galactic 
halo, and SWCX lines for modeling the background in the 
subsequent analysis, as well as the slope and normalization 
of the soft proton component. 


3 SPECTRAL ANALYSIS 

For the first time, we have enough X-ray photons to per¬ 
form spectral fitting in multiple annuli around an isolated 


spiral galaxy. We define nine concentric regions, whose lay¬ 
out is illustrated in Figure 4, along with an optical image 
of the galaxy. The hrst region is a circle of radius one ar- 
cminute centered at the nucleus of NGC 1961. The other 
eight regions are annuli of width one arcminute, centered 
at projected radii of I’.O, 1’.5, 2’.0, 2’.5, 3’.0, 3’.5, 4’.0, and 
4’.5 from the nucleus of the galaxy. Each region therefore 
overlaps with the adjacent regions. We chose this layout to 
improve the spatial resolution of our spectral analysis, and to 
test the robustness of our spectral modeling (i.e. overlapping 
regions should give consistent spectral results). However, the 
overlap means that adjoining spectra are not independent of 
one another, which we show graphically in subsequent plots 
through the X error bars spanning the radial range covered 
by each spectrum. We refer to these these regions as RO 
through R8, with the number increasing with projected ra¬ 
dius. RO covers the majority of the galaxy and its bulge, R1 
is dominated by the disk, and R2 by the outer disk. The 
other six regions capture the hot halo. 

We extract the spectra from each of these nine regions 
in each of our 30 observations. For each region, we perform a 
joint fit to these 30 spectra using the model described in sec¬ 
tion 2, with an additional set of components for the source 
emission. We use a redshifted photoabsorbed APEC -I- pow- 
erlaw to describe the source emission. We fix the redshift to 
the value for NGC 1961 {z = 0.013). We try to fix as many 
free parameters as possible before performing the spectral 
htting. We use the hts to the spectrum in the background 
region to get many of these constraints. The normalizations 
of the SWCX lines, the CXB, the Local Bubble component, 
and the Galactic halo are all fixed to the values obtained 
from the fit to the background region (scaled appropriately 
for each observation to account for the differences in angu¬ 
lar area). In our fiducial model, we also freeze the X-ray 
binary (XRB) component, the soft proton component, and 
the metallicity of the APEC component. We now discuss 
each of these choices in detail. 

The X-ray binary component is modeled as a power- 
law with redshifted photoabsorption from NGC 1961 as well 
as Galactic absorption. In our fiducial model we freeze the 
slope of this powerlaw at a reasonable value of F = 1.56 


(Irwin et al. 20031. For this model we also set constraints 


on the normalization of this powerlaw, by examining the 
surface brightness profile of the hard-band image (see next 
section). In the hard-band image, the hot halo of NGG 1961 
is expected to contribute negligible emission, so we can at¬ 
tribute all the observed emission to XRBs. We measure the 
background-subtracted hard-band flux in each region and 
use this flux (and an assumed F = 1.56 powerlaw spec¬ 
tral shape) to estimate the expected soft-band flux. This 
conversion is slightly uncertain, primarily due to uncertain 
background subtraction in the hard-band image and due to 
unknown absorption from NGG 1961 (since we do not con¬ 
strain the absorption from NGG 1961 until performing the 
spectral fitting), so we allow the normalization of the XRB 
component to vary slightly around the expected value in our 
fiducial model. 

In order to test that our modeling of the XRBs does not 
introduce any systematic effect, we also explore an alternate 
set of spectral models where we relax the above constraints. 
We let the normalization float as a free parameter. In region 
RO, we also let the slope vary (between F = 1.2 and F = 
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The Hot Halo Around NGC1961 5 



(a) Regions RO, R2, R4, R6, and R8 


(b) Regions Rl, R3, R5, and R7 


Figure 1. An optical DSS image of NGC 1961 (the same image is shown in both panels), along with our nine primary spectral extraction 
regions. Each region has a width of 1’ (16.6 kpc). Note that regions Rl, R3, R5, and R7 (shown in right) overlap with regions RO, R2, 
R4, R6, and R8 (left). The first region (RO) covers the majority of the galaxy and its bulge, while the outer disk dominates region R2. 
Regions R3-R8 capture the emission from the hot halo. 


2.0), and we tend to find that a harder slope than F = 
1.56 is preferred, probably due to a low-luminosity AGN. 
In regions Rl and R2, we let the slope vary between F = 
1.5 and F = 2.0 to allow for the presence of high-mass X- 
ray binaries (which are typically softer than low-mass X- 
ray binaries). In the other regions, where very little XRB 
emission is expected, we fix F = 1.56; in these regions the 
XRB component is typically unimportant, and the data are 
insufficient to constrain the slope. 

The soft proton component is also modeled with a pow- 
erlaw, though this powerlaw is not folded through the stan¬ 
dard instrumental response. We freeze the slope of this pow¬ 
erlaw to the slope measured in the background region (as 
[Snowden fc Kuntz||2014] note, it is generally reasonable to 
assume the spectral shape of the proton background is spa¬ 
tially constant across a given observation). In our fiducial 
model we also freeze the normalization of the soft proton 
component, using the ESAS proton_scale routine in order 
to convert the normalization from the measured value in the 
background region into the appropriate value for each re¬ 
gion for each observation. Again, in order to ensure that our 
constraints on this component do not systematically affect 
our conclusions, we also explore an alternate set of spectral 
models where we relax the above constraints. In this case, 
we relax them by allowing the normalization of the soft pro¬ 
ton background to vary by up to 100% in each observation, 
while keeping the slope fixed at the fiducial value. 

Finally, the metallicity is a key parameter of interest, 
but this is a very difficult quantity to measure observation- 
ally. We discuss this parameter in more detail below (in 
section 3.1), but in brief the difficulties arise due to a de¬ 
generacy between the metallicity and the normalization. In 


our fiducial model, we assume an intermediate value for the 
metallicity {Z = 0.5Zq), which reduces the magnitude of 
the possible error on the normalization (and the inferred 
density of the hot gas). To first order, our results can be 
scaled for a different assumed value for the metallicity by 
multiplying the density by the inverse square root of the 
fractional change in metallicity. However, again, we want to 
check that our assumption for the metallicity does not in¬ 
troduce a systematic error on the other components, so we 
explore an alternate set of spectral models with the abun¬ 
dance as a free parameter. 

We therefore have three sets of constraints we can toggle 
for our spectral modeling, yielding eight total spectral mod¬ 
els. The fiducial model uses all three of these constraints, and 
the other seven models relax one, two, or three of these con¬ 
straints. The dispersion between the results obtained from 
fitting with each different model gives us a way to estimate 
the systematic uncertainty in our results due to the use of 
simple spectral models for fitting to a complex astrophysical 
system. 

Finally, after finding the best-fit parameters for each 
model in XSPEC, we use the chain command to per¬ 
form a Markov Chain Monte Carlo (MCMC) search of the 
parameter space in order to account for degeneracies be¬ 
tween parameters and to make sure we properly sample the 
multidimensional space. XSPEC has two implementations 
of MCMC algorithms - Goodman-Weare and Metropolis- 
Hastings - and we explored both but due to the hard limits 
on the scaling factors for each instrument (allowed to range 
from 0.9-1.1) we achieved faster convergence and a higher ac- 
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ceptance fraction using the Metropolis-Hastings algorithnQ 
We initialized the chain using a diagonal covariance matrix 
with covariances taken from the XSPEC fits. We burn 10"^ el¬ 
ements before running the chain for 10® iterations, and apply 
a simple “simulated annealing” prescription for reducing the 
Metropolis-Hastings temperature of the fit as it progresses 
(we use an initial temperature of 5, and every 5000 steps 
we reduce the temperature by a factor of 0.9). This sam¬ 
ples the parameter space more fully before the chain begins 
to converge. We take the median value of each parameter 
(marginalized over the others) as the best-fit value and we 
also record the central 90% confidence interval of each pa¬ 
rameter in order to quantify the uncertainties. Some of these 
values are listed for the fiducial model in Table 2 below, for 
the regions where the emission measure of the hot gas com¬ 
ponent is at least 3cr above zero (note that the uncertainties 
listed in Table 2 bound the 90% central confidence interval, 
and are therefore larger than the Icr uncertainties). 

We also examine five larger annuli which are con¬ 
structed by combining the annuli shown in Figure 1. These 
annuli are 2’ in width, instead of the 1’ annuli shown in Fig¬ 
ure 1. These larger annuli have more photons and therefore 
yield better constraints, especially on the metallicity of the 
hot gas, but we find no systematic differences based on the 
different annuli sizes. 

In Figure 2 we present the resulting temperature profile, 
and in Figure 3 we present the metallicity profile. The black 
points (with 90% confidence regions) are the results for the 
fiducial model, and we plot the median results for each of the 
other seven models as well. The radii outside of which the 
normalization of the hot gas component falls below 3a are 
indicated as hatched shaded regions; in these outer regions 
the spectra fitting is not reliable. This occurs at r > 42 kpc 
for the 1’ annuli. For the 2’ annuli the hot halo is detected 
at >3cr in every annulus. 

The temperature in Figure 2 appears to decline slowly 
with radius. We can quantify the significance of this de¬ 
cline using the statistic, comparing a constant model to 
a linear model (so the two models differ by one degree of 
freedom). We fit to the spectral results for regions RO, R2, 
and R4, which are non-overlapping. We also include a sys¬ 
tematic error on each temperature, in quadrature, which 
is derived from the standard deviation of the eight best-fit 
values for each region from the eight different models we 
consider. Fven with the systematic error folded in, the lin¬ 
ear model is highly favored, with a Ax^ of 6.5 as compared 
to the flat model, giving 2.5cr evidence for a negative tem¬ 
perature gradient. As an alternative metric, since the error 
bars are asymmetric and we have the full probability dis¬ 
tributions from MCMC fitting, we also compute the Akaike 
information criterion for fiat and linear models. This gives a 
similar result, favoring the linear (negative) profile at 97.7% 
confidence (approximately 2.3a). On the other hand, for the 
metallicity, both methods show that a linear model is not 
at all favored over a flat model at the present level of sta- 


^ During our analysis, XSPEC was updated to introduce changes 
to the operation of the MCMC chain command. The MCMC re¬ 
sults in this work have been obtained using the newest version of 
XSPEC (v. 12.8.2q). 


tistical accuracy, and we conclude the metallically profile is 
consistent with being flat. 


3.1 On the low value of the metallicity 


The metallicity profile in Figure 3 shows strong statisti¬ 
cal evidence for sub-Solar metallicity throughout the hot 
gaseous halo. While the hot gas abundance is observed to 
be sub-Solar in some X-ray faint elliptical galaxies (e.g. Su 
|fc Irwin|2013[ ), metallicities as low as ours are still unusual, 
especially since a star-forming galaxy like NGC 1961 can be 
expected to have a higher supernova rate than a comparable 
elliptical. Such a low metallicity might suggest an external 
source (i.e. intergalactic medium) for the majority of the 
hot gas instead of an internal source. We think this result 
therefore warrants a bit more discussion. 

At temperatures around 0.6 keV, the metallicity is in¬ 
ferred from the ratio of the Fe L complex at around 0.7-0.9 
keV to the pseudocontinuum at around 0.4-0.55 keV. This 
procedure breaks down when the metallicity becomes high 
enough that line emission begins to dominate over the con¬ 
tinuum, so we performed simulations with XSPEC in order 
to verify that this is not a concern for these sub-Solar metal¬ 
licities, and to verify that the uncertainties returned by our 
MCMC modeling are of the expected order for the num¬ 
bers of photons in our spectra (approximately 10'* for the 1’ 
regions and twice as high for the 2’ regions). 

The reader can get a sense of the statistics in our spec¬ 
tra from Figure 4. For this figure, we have added the 20 
MOS spectra for each region to generate composite spectra. 
We emphasize that we do not fit models to these 
composite spectra; we fit to the individual spectra and 
propagate the backgrounds and angular areas separately for 
each spectrum. We show examples of the major model com¬ 
ponents in Figure 4 to illustrate their shape. The composite 
spectra show that the continuua are generally fairly well de¬ 
termined. The spectra behave roughly as expected as well, 
with the Fe L complex becoming increasingly weak as we 
move outwards in radius. Note also that we have not added 
the PN spectra to these composites; the PN spectra have 
roughly the same number of photons as the MOS spectra, 
improving our statistics by another factor of two. 

There are potential systematic errors, however. One 
possibility is the tendency of single-temperature fits to 
multi-temperature plasmas to yield anomalously low metal¬ 
licities ( Buote|2000 1, which we show in section 3.2 does not 
significantly affect our result. Another issue is the deter¬ 
mination of the continuum, and our model contains many 
components so it is important to check whether any of them 
may influence our measurement of the continuum. As we 
will show in the spatial analysis, the hot gas is the domi¬ 
nant component of the soft-band flux within the inner two 
or three arcminutes, after which the sky background and the 
QPB become the largest and second-largest components, re¬ 
spectively. The QPB is mostly expressed through the bright 
instrumental lines, none of which lie near the 0.4-0.6 keV 
continuum, so we do not expect it to affect the inferred 
metallicity. The sky background is modeled with a combi¬ 
nation of a power-law (for the CXB) and two APEC models 
(for the Calactic halo and the Local Bubble), each of which 
we fix based on the measurement in the background region, 
which begins 8’ from the center of the galaxy. A fluctuation 
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Table 2. Spectral Fitting Results 


region 

kT 

f neTifjdV 

Area 

/ d.o.f. 


(keV) 

( 10®2 cm-3) 

(arcmin^) 


RO 

n yq-|-0.02 
^'^-0.02 

96 3+1-1 

2.8 

1026.6/587 

R1 

Q 0O-I-O.O4 
O.DO_q Q3 

iO'^-0.8 

5.4 

646.7/471 

R2 

Q 0O-I-O.O8 

O.DO_o_o8 


7.7 

380.7/311 

R3 

0 

^•^-*--0.07 


9.2 

418.8/368 

R4 

0 30+0-25 

4 7+2.0 

10.1 

362.9/373 

R012 

rj ' 7 q-|- 0.02 
'*^-0.02 

34.4li;9 

10.5 

1122.8/819 

R123 

0 61+0-03 

21.71}! 

14.6 

1055.8/879 

R234 

^•of_o.io 

14.5!!! 

17.8 

956.6/877 

R345 

0 37+0-12 

10.7!}! 

21.7 

1012.8/955 

R456 

0 36+0-2® 

0-O0_Q QY 

5.9!!! 

26.3 

1138.8/1097 


Results for the hot gas component of the emission from NGC 1961, 
based on spectral fitting using our fiducial model. This model has the 
metallicity frozen at Z = O.SZq; for other choices of metallicity, the 
normalization should be scaled by a factor of approximately ^0.5/Z. 
The values and quoted errors are based on our MCMC chains. The 
best-fit values are the medians of the chain and the uncertainties 
bound the 90% central confidence region around the median. The 
area is determined from the BACKSCALE parameter in the spectrum, 
and shows the angular area of the fiducial spectrum for each region; 
this can be used to convert the emission measure into an average elec¬ 
tron density, as we do below. In each of the listed regions, the hot gas 
component of the spectrum is significant at more than 3cr. 


radius (kpc) 

10 20 30 40 50 60 70 80 



(a) 1’ annuli 


radius (kpc) 


10 20 30 40 50 60 70 80 



(b) 2’ annuli 


Figure 2. Temperature profile of the hot halo of NGC 1961, as measured using 1’ annuli (left) and 2’ annuli (right). Temperatures 
in each region are measured using eight different spectral models, as explained in Section 3. The best-fit values for the fiducial model, 
along with 90% confidence intervals, are shown in black. The best-fit values for the other seven models are shown in red and give a 
sense of the systematic uncertainty stemming from the choice of models. Confidence intervals for the red points are similar in size to 
the confidence intervals for the fiducial model, but are not shown for clarity. For points in the shaded hatched region (r > 42 kpc), the 
hot gas component in the model has less than 3cr significance; results in this region are not reliable. The temperature appears to decline 
slowly with radius. 


in these components on angular scales of several arcmin- 
utes would therefore lead to an incorrect assessment of the 
continuum in our source spectra. We can estimate the ex¬ 
pected magnitude of the CXB fluctuations using the results 
of Kolodzig et al. (in prep) for the XBootes field. In the 0.5- 
2.0 keV band on 4’ scales they find a typical size of CXB fluc¬ 
tuations of 5 X 10“^ ct^ s“^ deg“^ (Chandra ACIS-I counts). 
This corresponds to roughly 100 soft-band counts (XMM- 


EPIC MOS + PN instruments) in a 16 square-arcminute 
region observed for 100 ks, which is less than a percent of 
the counts in our spectra. It is possible that fluctuations be¬ 
tween 0.4 and 0.5 keV are more significant, however, since 
they include more Galactic halo and Local Bubble emis¬ 
sion, and the spatial variation of these components is not 
known. It is therefore not possible to firmly rule out the 
possible effect of background fluctuations from the CXB or 
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radius (kpc) 

10 20 30 40 50 60 70 80 



(a) 1 ’ annuli 


radius (kpc) 

10 20 30 40 50 60 70 80 



Figure 3. Metallicity profile of the hot halo of NGC 1961, as measured using 1’ annuli (left) and 2’ annuli (right). Temperatures in each 
region are measured using eight different spectral models, as explained in Section 3. The best-fit values for the fiducial model, along with 
90% confidence intervals, are shown in black. The best-fit values for the other seven models are shown in red and give a sense of the 
systematic uncertainty stemming from the choice of models. Confidence intervals for the red points are similar in size to the confidence 
intervals for the fiducial model, but are not shown for clarity. For points in the shaded hatched region {r ^ 42 kpc), the hot gas component 
in the model has less than 3(T significance; results in this region are not reliable. The metallicity is also essentially unconstrained in the 
1’ R2 annulus and in the outermost 2’ annulus. Overall the profile is difficult to measure precisely, but it is consistent with being fiat at 
a value around O. 2 Z 0 . 




Figure 4. Illustrations of stacked MOS spectra from various regions. The left plot shows stacked MOS spectra from each region, rebinned 
so that each point has a S/N of at least 10. Note that the hot gas component of the spectrum (roughly between 0.6 and 1 keV) becomes 
progressively less important at larger radii. In the right plot, we show the R2 spectrum as well as examples of representative spectral 
models for hot gas (red solid line), X-ray binaries (red dashed line), sky-|-instrumental background (thin black line), and the sum of these 
three components (thick black line). We stress that these stacked spectra and the models are for illustration purposes only. The actual 
spectral fitting is performed using simultaneous fitting to each of the 30 individual spectra for each region. 


Galactic backgrounds, although it would require significant 
fluctuations on these scales in order to bias our metallicity 
measurement. 

The SWCX and soft proton backgrounds also contribute 
to the 0.4-0.6 keV band. In the background region the 
SWCX component is about an order of magnitude below the 
CXB, however, and the SWCX background is not seen to ex¬ 
press significant spatial variations within an XMM-Newton 


field of view ( [Snowden & KuntzpOld] ), so we do not think 
this component is likely to affect the metallicity measure¬ 
ment. The soft proton background is also similarly low in 
the background region, and if we were incorrectly measuring 
it in such a way as to affect the inference of the metallicity, 
we would expect to see the inferred metallicity vary between 
the models where we freeze the soft protons and the models 
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where we allow some freedom in fitting this background; no 
such variation is observed between these models. 

An incorrect neutral Hydrogren column would also lead 
to incorrect estimation of the continuum. The Galactic Nh 
column was estimated from Dickey & Lockman (19901, but 


Kalberla et al. (2005 I also finds a similar value, indicating 


that the Galactic column is not particularly uncertain at 
this location. We checked the Planck CO maps as well and 
see no molecular gas in the direction of NGC 1961. 

Finally, NGC 1961 contains significant amounts of neu¬ 
tral gas in addition to the hot X-ray gas. It is therefore 
plausible that interactions between these phases would pro¬ 
duce charge exchange emission. Depending on the amount 
of this emission, and on the ratio of O vil to O vill emission, 
this effect might also change the ratio of the continuum to 
the Fe L complex. This sort of emission is known to be im¬ 
portant in the starburst galaxy M82 (Liu et al.|2011 Zhang 


|et al.|2014[ ) and in a few other cases | Li fc Waiig|20i3[ ), but 
a systematic study has yet to be performed so it is not yet 
possible to estimate its importance in NGC 1961. 

In summary, our inference of a low metallicity for the 
hot halo of NGC 1961 seems to be robust against basic 
systematic errors, but it is not definitive. We therefore in¬ 
clude a factor of Z/0.2Zq in subsequent figures to emphasize 
this uncertainty and show the effect of different choices on 
the proceeding conclusions. We also remind the reader that 
the Iron abundance in [Anders fc Grevesse| ( |1989| ) is about 
40-50% higher than in other abundance tables. Accounting 
for this difference makes our metallicity consistent with the 
0.2 —O. 3 .Z 0 metallicity which is ubiquitously observed in the 
outskirts of galaxy clusters (Bregman et al.|2010 Werner et 

g[2^. 


3.2 Two-temperature fits 

Here we explore the viability of two-temperature fits to the 
hot halo of NGC 1961, instead of relying on a single APEC 
model. We focus on the 2’ annuli in order to ensure sufficient 
photon statistics, and we use the models with frozen proton 
background and frozen XRB component. We initialize the 
two components at 0.3 and 0.8 keV and let the metallicity 
float. We therefore have three additional degrees of freedom 
(temperature, metallicity, and normalization of the second 
hot gas component) as compared to the one-T model, so 
for the second component to be statistically signihcant at 
p = 0.99, we require it to improve the by at least 11.35. 

The inner three regions (which overlap with the disk of 
the galaxy) show significant improvement with a 2-T model, 
while the outer two regions (which only cover the hot halo) 
do not show much improvement with a 2-T model. This 
seems intuitively reasonable, and one can imagine that one 
component could describe the hot halo and the other could 
describe the hot ISM of the disk. In general, the hrst com¬ 
ponent has a low metalllicity and contains most of the mass; 
its temperature, metallicity, and density behave like the 1-T 
hts shown in the previous section. This component can be 
associated with the hot halo of the galaxy. Note that, while 
the emission measures in Table 3 are larger than the fidu¬ 
cial emission measures, this difference stems from the lower 
metallicites in the 2-T fit as compared to the fiducial model 
(which has Z = 0.5Z©). The second component is more 
poorly constrained, but generally has a higher metallicity 


and a lower emission measure. This component can be asso¬ 
ciated with the hot interstellar medium of the galaxy. The 
2-T model therefore seems to broadly support our picture of 
the hot gas in and around NGC 1961. Unfortunately there 
are not sufficient photons to do a more detailed analysis of 
the ISM of this galaxy at present. 


4 SPATIAL ANALYSIS 

In this section we examine the surface brightness profile of 
the X-ray emission around NGC 1961. We analyze the soft- 
band and hard-band images generated in Section 2, as well 
as model images of the various background components. The 
ESAS mos_back and pn_back routines generate model par¬ 
ticle background images automatically, and we use the ht 
parameters from our ht to the background at large radii 
as inputs to the the ESAS soft_proton and swcx routines 
in order to generate models of the soft proton background 
and the SWCX backgrounds for each observation and each 
detector (these two backgrounds are only generated for the 
soft band, however). For each type of image, we use the 
merge_comp_xmm routine to combine all 30 images together 
into a single merged image. We also use this routine to com¬ 
bine the 30 exposure maps together to generate a merged 
exposure map. We present the merged, exposure-corrected, 
soft-band image in Figure 5, which was created with the 
merge_comp_xmm routine and has the QPB, soft proton, and 
SWCX backgrounds all subtracted. Note that this rou¬ 
tine applies no weighting between the MOSl, MOS2, 
and PN detectors, but instead adds the images (in 
counts) directly; this has the effect of weighting the PN 
image more than the MOS images due to the larger effective 
area of the PN detector. 

We divide each image by the merged exposure map and 
construct surface brightness profiles in Figure 6. The net 
soft-band signal flattens to a constant value at large radii; 
we use this mean value as the value for the X-ray background 
in this field over the soft band. The X-ray background (the 
combination of the CXB, Galactic halo, and Local Bub¬ 
ble) is the largest background component, followed by the 
particle background, the soft proton background, and h- 
nally the SWCX background. The sky background estimated 
from Figure 6 is consistent with the 3/4 keV background 
for the same region measured from ROSAT as well. Note 
that the particle background and the proton background 
appear to turn upwards at large radii; this is because these 
backgrounds are not focused by the telescope’s mirrors and 
therefore do not suffer the same vignetting as the X-ray 
backgrounds. Dividing by the exposure map therefore over¬ 
estimates these backgrounds at large radii. This should not 
be a major issue for our analysis, since it does not begin to 
matter until radii of 100 kpc or more, which is outside the 
range within which we can measure the hot halo. 

In these profiles we also show, but do not subtract, the 
estimated contribution of X-ray binaries to the soft band. 
We estimate this contribution in two different ways. One 
method (the cyan line) uses the hard-band X-ray image and 
model backgrounds. We determine the surface brightness 
profile of the hard-band emission and multiply this profile 
by a factor of 1.67 (the scaling into the soft-band summed 
MOS1-I-MOS2-I-PN counts for a F = 1.56 powerlaw with 
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Table 3. 2-Temperature spectral fits 


region 

Ti 

Zi 

f neTifjdV 

72 

Z 2 

f UenjjdV 

Ax® 


(keV) 

(Ze) 

(10®2 cm-3) 

(keV) 

(Zq) 

(10®2 cm-®) 


R012 

n 75+0 09 
' O_0.04 

0 19+0-07 

71 e + 20.4 
'-*-•0-24.8 

'-'•^0-0.07 

1 54+0-S5 

-L-^^_0.82 

6.7l« 

30.8 

R123 

0 74+0-10 

n 1 Q+O-io 
u.io_o 08 

36 . 81 ?®:® 

n 90+0-07 

^•^o_0.09 

-L-yD_1.64 

er 0+21.3 

0-'^-2.6 

33.3 

R234 

Q 05 + 0.12 

'J-'50_o.07 

^•^O_o.i8 

35.4l«;® 

J--'0'0_0.35 

-| 70+2.67 
''^-1.18 

l.S+a 

24.3 

R345 

n “^s+O-1^ 
lJ-'3O_0.05 

0.27t°it 

1 7 q+18.0 

0 02"'"®'^^ 

2 . 301 ?;?? 

1.2ll?i X 10® 

2.2 

R456 

'J-'50_o.08 



0 02"'"®'^^ 

o.oilHl 

5.5l|;? X 10® 

0.2 


Two-temperature fits to the 2’ regions. The final column shows the improvement in the goodness of 
fit parameter for the 2-T model relative to the 1-T model; a value of at least 11.35 corresponds to an 
improvement significant at more than 99% confidence. The 2-T model is statistically favored in the inner 
three regions, which include the disk of the galaxy, but offers no significant improvement in the outer 
regions. 



Figure 5. Adaptively smoothed merged XMM-EPIC image of the 0.4-1.25 keV emission from the region around NGC 1961. This 
image has been exposure-corrected and the estimated particle background, soft proton background, and SWCX background have been 
subtracted; point sources have also been masked. The CXB, Galactic halo, and Local Bubble have not been subtracted and these 
backgrounds produce the uniform background which fills the field. The white arrow indicates 3’, or about 50 kpc at the distance of NGC 
1961, and the black cross denotes the center of the optical emission from the galaxy. The hot halo is visible by eye to several arcminutes 
and can be studied at larger radii through surface brightness profiles. 


Galactic absorption) to generate the estimated XRB profile 
in the soft band. For the other method (magenta line), we 


use the K-band image of this galaxy from 2MASS (Skrut- 
skie et ar]|2006 1 to estimate the K-band surface brightness 


profile, then convert to the soft band using the scaling rela- 

poTlI (this 


tion for LMXBs from Boroson et al. 


is also described in more detail in Anderson et al. 20131. 
These two methods agree reasonably well, both showing a 
subdominant but somewhat uncertain XRB contribution in 


the central region, which falls off quickly beyond about 1.5 
arcmin (23 kpc) and becomes insignificant at larger radii. 

We also verify that our results are consistent between 
the spectral and spatial analyses. This is important to check 
(Anderson and Bregman 2014), especially with the complex¬ 
ity of our background model. To check this, we plot the sur¬ 
face brightness profile of the hot gas as determined from each 
technique, in Figure 7. The black line shows the result of 
the spatial analysis, which is the remaining soft-band emis- 
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radius (kpc) 

2 5 10 20 50 100 200 



Figure 6 . Background-subtracted surface brightness profile of the 0.4-1.25 keV emission around NGC 1961. The black points are the 
background-subtracted data, showing the remaining emission after subtracting: (red) particle background, (green) soft proton background, 
and (blue) Solar wind charge exchange background. The upper set of black points include the contribution of sky background, which is 
assumed to be flat across the field and fit with the dashed black line. The lower set of black points has the estimated sky background 
removed as well, leaving only the emission we attribute to NGC 1961. Finally, the dotted cyan and magenta lines show the estimated 
contribution of X-ray binaries in NGC 1961, using two different methods. These are not subtracted from the profile but are clearly 
subdominant in the soft band. 


sion after subtracting the estimated CXB, quiescent particle 
background, soft proton background, Solar wind charge ex¬ 
change background, and XRB emission as estimated from 
the hard-band image (the cyan line in Figure 6). We restrict 
this plot to the MOS instruments, since the PN instrument 
has a different effective area. The red points are the results 
of the spectral analysis (using the fiducial model), where 
we have converted the APEC component describing the hot 
halo (including Galactic and intrinsic absorption) into a soft- 
band count rate for the MOS detectors. 


The agreement is excellent within the ~ 42 kpc where 
the spectral fits are robust. In the outer spectral regions, 
where the significance of the hot gas component is less than 
3 (t, the spectral model predicts far lower surface brightness 
than the observed spatial profile. However, here the hot gas 
emission comprises less than about 10% of the total soft- 
band signal, and it is not possible to say with certainty what 
the true emission looks like at such low surface brightnesses. 
We discard these outer regions from the joint analysis below. 


5 PHYSICAL PROPERTIES OF THE HOT 
HALO 


Now that we have performed both a spectral analysis and a 
spatial analysis, and shown that they give consistent results, 
in this section we explore the physical properties of the hot 
halo. First, we deproject the spectral results in order to es¬ 
timate the hot gas density profile. Conceptually, we follow a 
similar procedure to that of |Churazov et ah (20031, although 
these data have much lower S/N than their observations of 
the Perseus cluster so we make a few modifications to that 
procedure. 


First, we generate a simple estimate for the emission 
measure of the hot halo at very large radii. To do this, we 
fit a power-law to the hot gas component of the observed 
surface brightness profile at large radii, finding a logarith¬ 
mic slope of 3.5 for the surface brightness as a function of 
radius. This corresponds to a slope of 2.25 for the density as 
a function of radius, which is equivalent to /3 of 0.75 in the 
standard /3-model. This conversion between surface bright¬ 
ness profile and density profile assumes that the emissiv- 
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radius (kpc) 



Figure 7. Spatial (black points) and spectral (red points) com¬ 
parisons for the MOS soft-band surface brightness profile at¬ 
tributed to the hot gaseous halo. For comparison, we also show the 
soft-band surface brightness profile of the full image (i.e. instru¬ 
mental + sky backgrounds; dashed line) and the surface bright¬ 
ness profile from a model hot halo with a uniform density of 
200pcf2f) (dotted line). Black points with open circles denote neg¬ 
ative values which have been multiplied by —1 to allow plotting 
with a logarithmic y-axis. Red points with open circles are spec¬ 
tral fits where the significance of the hot gas component is below 
3 (7. The results from the spectral and spatial methods largely 
agree within r £ 42 kpc. At larger radii, where the surface bright¬ 
ness of the hot gas is below about 10% of the total signal, we 
cannot tell which is correct: the profile from the spatial method, 
the profile from the spectral method, or neither. 


ity does not vary with radius. In our 0.4-1.25 keV band, 
the emissivity does not change with temperature by more 
than 10% for temperatures between 0.3 keV and 1 keV. If 
the metallicity varies, the emissivity can also change, but 
for this galaxy the projected metallicity profile is consistent 
with being fiat. We use a fiducial temperature of 0.4 keV 
and a fiducial metallicity of 0.2 Zq to convert from surface 
brightness into density. 

Next, we convert the results of the spectral fitting into 
estimates of the emission measure in each region. This is 
trivially derived from the normalization of the APEC com¬ 
ponent of the spectrum, using our assumed distance of 58.0 
Mpc to NGC 1961. We rescale the emission measure in each 
region assuming the fiducial metallicity of 0 . 2 ^ 0 . While we 
adopted a metallicity of 0.5Zq when performing our spectral 
fitting, the results of the spectral fits with floating metallic¬ 
ity suggest that a value of O. 2.^0 is more reasonable for this 
galaxy. We therefore adopt a fiducial metallicity of 
O. 2 Z 0 for the remainder of the analysis. It is simple to 
scale our results into a different metallicity, however: at the 
level of accuracy we can measure for this galaxy, a good ap¬ 
proximation is that the emission measure is inversely propor¬ 
tional to the value of the metallicity, and the electron density 
is proportional to the inverse square root of the metallicity. 

We treat the spectral fits to the 1’ bins and the 2’ bins 
separately, and since the regions overlap we also separate 
each set into two groups. We therefore have four indepen¬ 
dently determined profiles, based on regions R0-R2-R4, Rl- 
R3, R012-R345, and R123-R456 respectively. As noted in 
section 4, we discard the R5-R8 regions since the normal¬ 


ization is so poorly constrained from the spectral fitting in 
these regions. 

Finally, for each annulus, we subtract the expected 
emission measure from exterior annuli (EMext). We divide 
the remaining emission measure by the volume of the annu¬ 
lus, scaled by the fraction of the field of the annulus covered 
by the fiducial spectrum (using the area A listed in Table 
2). We also include a factor of 0.83 to convert from uh into 
rie assuming a standard Helium abundance. Putting this all 
together, and using a distance of 58.0 Mpc and angles in 
units of arcminutes, the expression for the average electron 
density is (see also McLaughlin] 1999 who derives a more 
general form of this equation): 


Ue = 8.00 X 10 ^ cm ^ X 


f UeUndV — EMext 62 — of 
10 ®® cm“® 


A (01-Of) 

( 1 ) 

These values are listed in the final column of Table 2, 
and displayed graphically in Figure 8 . The systematic errors 
seem to be larger than the Icr statistical uncertainties, but 
together they are still only at about the 10% level. The un¬ 
certain metallicity is by far the largest source of uncertainty. 
We also note that in the inner three regions, deviations from 
hydrostatic equilibrium may be expected due to the presence 
of the galactic disk. There is some evidence for this in the 
preference for 2 -temperature fits to these spectra, but the 
underlying hot halo component remains dominant in these 
regions, and there is no evidence of the disk in the X-ray 
image (Figure 5). 

In Figure 8 we also plot the best-fit hot halo electron 
density profile for this galaxy, as measured by B13. They 
used the modified /3-model profile of Vikhlinin et al. (2006 1 
to parameterize the surface brightness profile, and assumed 
a constant metallicity of 0 . 12^0 (relative to the abundance 
tables of Grevesse fc Sauval|1998 1 . We also compute a “cor¬ 
rected” version of their density profile, which is multiplied 
by a factor of 0.64 to account for the different metallicity 
and the different abundance table relative to our analysis. 
Overall the agreement is very good between our deprojected 
profile and their corrected best-fit profile. The behavior of 
the profile at larger radii is extremely important, however, 
and it is not clear whether their parameterization can be ex¬ 
tended to larger radii. Improved observations are necessary 
in order to understand the behavior of the hot gaseous halo 
within a larger fraction of the virial radius. 


5.1 Pressure and Mass Profiles 

We can estimate an electron pressure profile for the hot gas, 
which is the product of the electron density and the temper¬ 
ature. Unfortunately, we do not have a deprojected tempera¬ 
ture profile, and the hot gas component of the spectra is sub¬ 
dominant beyond a projected radius of 2 arcminutes, so we 
cannot get robust results by subtracting scaled spectra from 
one another and fitting the remainders, as one can do for de¬ 
projection in the high S/N regime. However, our projected 
temperature profiles do not show significant gradients, and 
we can produce an approximate estimate for the reprojected 
pressure profile by multiplying the reprojected density pro¬ 
file by the projected temperature profile. We propagate the 
uncertainties on the temperature and the density into the 
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(a) 1’ annuli (b) 2’ annuli 

Figure 8. Deprojected electron density profiles of the hot halo of NGC 1961. Since the data points overlap, each plot shows two 
independent profiles, derived from non-overlapping points. The two independent profiles are perfectly consistent with one another. As 
in the other figures, the black points show the results for our fiducial model and the red points show the results for the other spectral 
models, and the error bars are Icr around the median for the fiducial model. For comparison, the density profile for this galaxy measured 
by B13 is displayed as well, along with a “corrected” version of their density profile rescaled to match our fiducial metallicity of 0.2Zq. 
The corrected profile matches our results well. 



Figure 9. Approximate electron pressure for the hot halo of NGC 
1961, based on the temperature profile and the deprojected elec¬ 
tron density profile (see section 5.1). Errors are Itr. 


total uncertainty, and neglect the additional uncertainty in¬ 
troduced by using a projected temperature profile instead 
of a deprojected temperature profile (this should be smaller 
than the statistical uncertainties on the temperature, how¬ 
ever). The resulting profile is shown in Figure 9. 

Now that we have deprojected electron density and elec¬ 
tron pressure profiles, we can estimate the total mass profile 
of the galaxy. This is derived from the temperature, elec¬ 
tron density, and the gradient of the total pressure at each 
radius. We assume jr = 0.61 and we assume the total pres¬ 
sure is 1.91 times the electron pressure. This factor is appro¬ 
priate for Solar metallicity gas with the [Anders fc Grevesse| 
(19891 abundances, but does not vary by more than 1% for 


tween the points in our pressure and density profiles in or¬ 
der to calculate the gradient, and extrapolate to larger and 
smaller radii based on the distance-weighted slope of the 
nearest data points (for more details, see Appendix B of 


Churazov et al. 20081. The effective circular velocity cor¬ 


responding to this derived mass profile is shown in Figure 

10 . 

For comparison, we also plot a number of other con¬ 
straints on the circular velocity profile of this galaxy. |Haan| 


et al. (2008 \ measured the HI circular velocity out to 43 


kpc, with independent measurements for the receding and 
approaching sides of the disk. They measure an inclination 
angle of i = 42.°6 ± 4°.0 (close to the 47° listed in Hyper- 
Leda) which they use to correct their measurements. The 
approaching side of the disk decreases towards zero veloc¬ 
ity near the center, but the receding side of the disk shows 
a roughly fiat curve towards the center. At smaller radii, 
CO 1-0 measurements are also available from ICombes et al.l 
(20091, which are more consistent with the behavior of the 


receding side of the disk when corrected using i = 42.°6. 


other reasonable abundances. We neglect magnetic support 
or other forms of non-thermal pressure. We interpolate be- 


Combes et al. (20091 make an argument for a higher 
inclination angle, however, using i = 65° as a rough value 
which we adopt here for comparison. This value would re¬ 
duce the inclination correction and the inferred total mass 
for the galaxy. Applying this inclination correction seems to 
introduce some tension between the cold gas measurements 
and our X-ray measurements, however. 

We consider a third mass estimate using the K-band im¬ 
age to derive the stellar contribution to the potential. This 
yields a lower limit to the total circular velocity. We extract 
the stellar profile from the K-band image using ellipses for 
both i = 42.°6 and i = 65°, but the results are nearly iden¬ 
tical in both cases so we only plot the former for simplicity. 

Using the assumed M/L ratio of 0.6 from Section 1, 
which is roughly the expectation based on [Bell fc de Jong| 
(20011 and Chabrier (20031 for this galaxy, we see that the 
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Figure 10. Rotation curves for NGC 1961. Data points are in¬ 
clination corrected assuming i = 42.°6 or i = 65°, as indicated. 
HI data are form |Haan et aL] | |2008| l for the approaching (circles) 
and receding (squares) sides of the disk. CO 1-0 data are from 
[Combes et al.| ( |2009^ . The red line is our estimate of the effective 
circular velocity based on our X-ray data, assuming hydrostatic 
equilibrium. The line is shaded thick over the region where the 
rotation curve is constrained by data, and dotted where the curve 
is extrapolated. The blue lines show the approximate contribution 
of the stars in NCC 1961, based on the K-band image assuming 
either a Chabrier (solid line) or Salpeter (dashed line) IMF. In 
general the i = 42.°6 model seems to be preferred by the X-ray 
data, and the rotation curve seems to be roughly flat within at 
least 42 kpc. 


stars contribute the majority of the mass in the central re¬ 
gion of the galaxy, and are also in tension with the gas- 
dynamical measurements corrected assuming i = 65°. If we 
instead consider the Salpeter (19551 M/L ratio, the tension 
increases further. The Haan et al.| ( |20()8 ( inclination angle 
gives results which seem consistent among gas-dynamical 
measurements, X-ray measurements, and the stars. In the 
regimes where each curve is reliable, the agreement between 
the rotation curves is good, and points to a roughly flat 
rotation curve within at least 42 kpc. 

Finally, in Figure 11 we plot the enclosed stellar mass, 
hot gas mass, and total mass for NGC 1961. We use the 
aforementioned M/L ratio of 0.6 for the stars, so that the 
total stellar mass is about 3 x 10^The rotation curves 
only extend to about 42 kpc, so we extrapolate the to¬ 
tal mass profile to larger radii. We assume that the dark 


matter approaches an NFW profile at larger radii (Navarro 


et al. 19971, such that the virial mass of this galaxy is 
1.3 X IO^^Mq, extending to a virial radius of about 490 kpc. 
This agrees well with the estimate of B13, who estimated 
a virial radius of 470 kpc for this galaxy based on compar¬ 
ison with cosmological simulations. Still, it is a very crude 
estimate of the total mass, and we will therefore not use the 
extrapolated mass for any precise calculations. 

Within 10 kpc, the stellar component seems to be dom¬ 
inant, but at larger radii the system quickly becomes domi¬ 
nated by dark matter. Within 90 kpc, the hot halo contains 
less than a tenth of the mass in the stars. At this radius, 
the baryon fraction (mass in stars -|- neutral Hydrogen -|- 
hot halo) is close to the Cosmic fraction. Extrapolating to 
larger radii, the hot halo may become comparable to the 


Figure 11. Enclosed mass profiles for NCC 1961. The green line 
shows the total mass, as inferred from the X-ray observations us¬ 
ing the assumption of hydrostatic equilibrium (section 5.1). The 
dashed extension at large radii is an NFW profile fit to the ob¬ 
served portion of the mass profile, as described in the text. The 
blue line is the stellar mass, as inferred from the 2MASS K-band 
image of this galaxy using a M/L ratio of 0.6. The red line is 
the mass of the hot gaseous halo, including la uncertainties as 
the shaded red region. Each line is styled in boldface over the 
regime where the mass is measured, and dotted where the profile 
is extrapolated. 


mass of the stars, but the sum of these components is less 
than a third of the expected baryon content for the system. 
This is an illustration of the problem of missing baryons 
from galaxies (see section 6.1). 

Our measured electron pressure prohle is also useful for 
predicting the thermal SZ signal from this galaxy. The ther¬ 
mal SZ effect is proportional to the volume integral of the 
electron pressure, which we express using the Compton y- 
parameter. The Compton y parameter for our pressure pro¬ 
file, integrated over a sphere with radius 42 kpc, is 2 x 10“® 
arcmin^. The typical radius for comparison with other stud¬ 
ies is TZsoo, which for our profile is ~325 kpc (19.6’) for 
NCC 1961 (and the corresponding M500 is 1 x lO^^M©). If 
we extrapolate our pressure profile out to this radius, the 
integrated y parameter increases to 1 x 10“® arcmin^. These 
values are below the sensitivity limits of Planck SZ catalogs 


for an object of this angular size (Planck Collaboration et 
|al.||2015^ |Khatri||2015| ) and the projected location of NGC 
1961 is fairly close to the Galactic disk, which is a more 
difficult location in the sky to isolate the SZ signal. 

We can also compare this value to the stacking mea¬ 
surements of Planck Collaboration et al. (20131. For galaxies 
with log M * /Mq = 11.5, Planck Collaboration et al. (20131 
estimate Fsoo ~ 4 x 10“® arcmin'^, where the tilde indicates 
that the value has been normalized to an angular distance of 
500 Mpc. Using the distance of 58 Mpc appropriate for NGC 
1961 yields Fsoo « 3 x 10“^ arcmin^. This is about 300 times 
larger than our estimated Isoo for this galaxy. Based on the 
Planck selection function data, for a galaxy with the angu¬ 
lar size of NGC 1961, the Planck map should be sensitive 
to such a signal at 5 (j or higher with around 90% complete¬ 
ness, so the non-detection of this galaxy suggests it does not 
have the mean thermal SZ signal for a galaxy of its stellar 
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Figure 12. Approximate entropy profile for the hot halo of NGC 
1961, based on the temperature profile and the deprojected elec¬ 
tron density profile (see section 5.2). The entropy rises with ra¬ 
dius, as expected for a dynamically stable halo. The central en¬ 
tropy is very close to the transition between cool-core and non- 
cool-core systems. 


mass. On the other hand, based on the halo mass of NGC 
1961, Planck Collaboration et al. (20131 predicts Ysoo of just 
1 X 10~'‘ arcmin^, which is undetectable in the Planck map. 
While this is only one galaxy, so it is difficult to draw strong 
conclusions, the implication is that halo mass is a better 
predictor of the SZ signal than the stellar mass, a point to 
which we return in section 6.2. However, the predicted Y 500 
for a galaxy with the halo mass of NCG 1961 from |Planck| 
|Collaboration et al.| ( |2013[ ) is still an order of magnitude 
above our predicted Y 500 value (which admittedly involves 
an extrapolation of nearly an order of magnitude in radius 
beyond our outermost spectral bin). More work needs to be 
done in order to understand this discrepancy. 


5.2 Entropy Profile and HI Emission 


We can also construct the entropy profile for the hot halo of 
NGC 1961. This is given hy K = , and propagates 

uncertainties on both T and ue just like the pressure profile 
(Figure 9). The entropy profile is shown in Figure 12. 

Within the uncertainties, the entropy profile is positive 
in agreement with the Schwarzschild stability criterion. The 
central value is about 30 keV cm^, which is the transition 


value noted by Cavagnolo et al. (20081 for brightest clus¬ 


ter galaxies (BCGs). BCGs below this value typically have 
bright Ha emission and star formation, and generally cor¬ 
respond to cool-core systems ( Cavagnolo et al^|2009 l, while 
higher entropy values indicate that the gas is able to strat¬ 


ify before cooling instabilities grow (McCourt et al. 2012 


Sharma et al. 2012 Caspar! et al. 2012 Voit et al. 20151. 


The low entropy in the core of NGC 1961 therefore leads 
one to expect similar behavior for this galaxy. Indeed, 21 
cm HI observations of this galaxy ( Haan et al.||2008 l show 
widespread, patchy HI emission around this galaxy, extend¬ 
ing tens of kpc to the NW. This extended HI emission has 
a low velocity dispersion (of order 20 km/s) and appears to 
be co-rotating with, and accreting onto, the disk (see Figure 
6 of Haan et al.|2008 1. This extended HI appears to contain 


about 10% of the total neutral Hydrogen content of the sys¬ 
tem, or about 5 x 10 ®Mq of neutral gas. In order to check if 
this neutral gas can be explained by condensation from the 
hot halo, we compute the cooling time and the dynamical 
time as a function of radius. The cooling time is given by: 

3/2nfcT _ 3/2fcT x 1.91 
Ane{n — rie) Arie x 0.91 

where n is the total particle density, and we assume n = 
l.Olrie as in section 5.1. The dynamical time is given by 


3tt 

16Gp 




AGM{< R) 


( 3 ) 


At every radius where we have data, the cooling time is 
about two orders of magnitude higher than the dynamical 
time, which is an order of magnitude away from the thresh¬ 
old value predicted from simulations for cooling instabilities 
to become important. At 40 kpc, for example, the cooling 
time is 13 Gyr and the dynamical time is 160 Myr. There¬ 
fore it seems unlikely that condensation from the hot halo 
is the source of most of the neutral Hydrogen seen around 
NGC 1961. 


Another possibility, proposed by Shostak et al. (19821 


is that the HI is the result of ram pressure stripping. NGC 
1961 lies several hundred kpc to the NW of a poor group 
of galaxies, and the HI appears extended in the direction 
away from the center of this group, so it appeared plausible 
that NGC 1961 could be falling into this group and losing 
its ISM through interactions with a hot intragroup medium. 
However, the X-ray halo of NGC 1961 shows no signs of such 
interactions, and there is also no X-ray emission associated 
with a larger intragroup medium in Figure 5. A ROSAT 


analysis by Pence & Rots (19971 also found no sign of an 


intragroup medium emitting in the X-ray band. 

We think the extended HI feature is more likely to be a 
sign of accretion from the intergalactic medium. This would 
explain its filamentary appearance, and it seems logical that 
an unusually massive spiral galaxy would require an unusu¬ 
ally large supply of intergalactic gas. Moreover, the velocity 
field of the HI stream appears identical to the velocity field 
of the HI disk (see Haan et al. 2008 Figure 6), which was also 


pointed out by Shostak et al. (19821. Alignment of the angu¬ 


lar momentum of the accreting gas and the galactic disk is 
the natural expectation for accretion from the IGM (Mo et 
al.|1998 1. Such alignment is less likely for a minor merger, so 
we think accretion is a better explanation than a merger for 
this galaxy, but a merger or accretion of a satellite galaxy is 
also a possibility. 

It is interesting that, although the entropy profile looks 
like a classic cool-core system, the temperature of the hot 
gas actually increases in the center by factor of 50%. While 
cool cores are typical in galaxy groups and clusters, negative 


gradients are common for galaxies (Humphrey et al. 2006 


Diehl & Statler 20081, although X-ray faint ellipticals can 


show flat or positive gradients as well (Fukazawa et al. 20061. 
NGC 1961 is an isolated galaxy with the mass of a galaxy 
group, but it seems that the halo mass is not the determinant 
of the shape of the temperature profile for this system. We 
speculate that feedback from the galaxy itself (either stellar 
feedback or a previous episode of AGN feedback) may be 
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responsible for heating the gas within the central few tens 
of kpc in this system. 


6 DISCUSSION 


6.1 Missing Baryons from NGC 1961 

In this section we return to one of the questions discussed in 
the introduction: where are the missing baryons from NGC 
1961? The deep observations in this study allow us to con¬ 
strain the density profile of the hot halo out to about 42 
kpc, which is about 9% of the virial radius, and it is notable 
that we see no evidence for any “flattening” in the density 
profile over this range. 

It is still difficult to rule out a putative fiat component 
at larger radii, however. As Figure 7 shows, the spectral 
and spatial techniques begin to disagree at radii larger than 
about 50 kpc, where the hot gas surface brightness is be¬ 
low about 10% of the total soft-band surface brightness. We 


are unlikely to improve upon this limit (see, e.g. Vazza et 
al. 2011| Reiprich et al. 20131 with the current generation 
of instrumentation. In Figure 7 we also show the predicted 
surface brightness from a uniform density cloud of gas with 
kT = 0.5 keV, Z = 0.2Zq, and n = 200pcBi) extending 
out to the virial radius. This is the flattest possible density 
profile that contains the entire baryonic content of the sys¬ 
tem. Unfortunately this limiting case falls below our surface 
brightness sensitivity, and so we cannot confirm or rule out 
the possibility of some sort of extremely flat halo becoming 
important at larger radii. 

It is also important to emphasize the uncertainties in 
this discussion. We can detect the hot halo to nearly 20% 
of the virial radius from our surface brightness analysis, and 
measure the density to about 10% of the virial radius with 
good accuracy, but the vast majority of the virial volume 
is unconstrained and depends on power-law extrapolations. 
Moreover, the total mass of NGC 1961 is also not well- 
constrained, since it also depends on similar extrapolations 
to the observed rotation curves. Finally, since we find the 
stars have a similar mass as the hot halo for this galaxy, the 
M/L ratio of the stars is also somewhat important. Our es¬ 
timates of each of these quantities are unlikely to be wrong 
by a factor of 2 or more, so the problem of missing baryons 
remains for NGC 1961, but the exact baryon fraction should 
be considered extremely uncertain. 

Our measurements of the pressure in the hot halo can 
also be useful for understanding the baryon budget of NGC 
1961. Observations with the COS instrument on HST show 
ultraviolet absorption around L* galaxies from many dif¬ 
ferent elements, generally with a high covering fraction 
(Tumlinson et al. 2011 Stocke et al. 2013| Werk et al. 
2014). Photoionization modeling suggests that this absorp¬ 


tion comes from ionized clouds, with characteristic temper- 

-1-3 „„-3 


and a volume filling fraction around 10% ( Werk et al.|2014 |. 
Such clouds cannot exist in pressure equilibrium with the 
hot halo around NGC 1961, but at larger radii (several hun¬ 
dred kpc) if the hot halo pressure continues to decline these 
clouds would be in equilibrium. If they have densities at the 
high end of their estimates and a 10% volume filling fraction 
beyond 150 kpc, such clouds could contain essentially all of 



Lk/L^ (t (km/s) 


Figure 13. 0.5-2.0 keV luminosity of the hot gas component 
of nearby galaxies, plotted against (left) K-band luminosity and 
(right) ID stellar velocity dispersion. Blue and black points 
correspond to field elliptical galaxies taken from the literature 
( [Mulchaey &: Jeltema||2010[ [Boroson et al.||2011[ |, and the black 
shaded region is the best-fit relation with Icr uncertainties from 
[Mulchaey fc Jeltema| l |2010[ |. The red star is NGC 1961, using 
measurements from this work. We have converted the peak circu¬ 
lar velocity from the HI profile into an effective velocity dispersion 
by dividing by \/3. 


the missing baryons from the galaxy. However, such clouds 
are unlikely to condense from the hot halo, given the high 
ratio of the cooling time to the dynamical time. If they orig¬ 
inate from the galaxy, they would therefore have to survive 
the passage through the inner 100 kpc or more of hot gas. 


6.2 Comparison with Other Galaxies 


We measure the X-ray luminosity of the hot halo for this 
galaxy from the surface brightness profile (Figure 6). We 
convert from 0.4-1.25 keV counts into 0.5-2.0 keV flux as¬ 
suming an APEG model with a fiducial temperature of 
0.5 keV and a metallicity of 0.2 Zq, although for reason¬ 
able values of Z and T the conversion factor is constant to 
within 5%. The hot gas luminosity within 50 projected kpc 
is 8.9 ± 1.2 X 10'*'^ erg s“^. Integrating out further yields a 
total hot gas luminosity of 1.5 ± 0.4 x 10^^ erg s“^. This is 
slightly higher, but given the much larger uncertainties the 
difference is not very significant. 

In Figure 13 we compare the 0.5-2.0 keV luminosity of 
NGC 1961 (within 50 kpc) to various scaling relations for 
early-type galaxies. The left panel shows this galaxy on a 
plot of Lx versus Lk, comparing against isolated galaxies 


from Mulchaey & Jeltema (20101 and from Boroson et al. 


(20111. NGC 1961 lies below their best-fit relation by about 


an order of magnitude. The scatter in these relations is also 
nearly an order of magnitude, so we caution against drawing 
strong conclusions from a single object, but this is consistent 
with the general observation that the X-ray halos of spiral 
galaxies appear less luminous than elliptical galaxies. 

NGC 1961 is much more consistent with relations be¬ 
tween Lx and the total matter content. In the right panel 
of Figure 13 we show the comparison between Lx and the 
velocity dispersion of early-type galaxies from |Boroson et al.| 
(2011). Since NGC 1961 is rotationally supported, we place 


this galaxy in the Figure by dividing the peak inclination- 


© 0000 RAS, MNRAS 000, 000-000 
















































The Hot Halo Around NGC1961 17 




radius (kpc) radius (kpc) 


Figure 14. Temperature, abundance, density, enclosed gas mass, 
entropy, and enclosed total mass profiles for NGC 7796 (green), 
(blue), and NGC 4455 (red), including Icr error regions, from 
O’Sullivan et al. (2004) and O’Sullivan et al. (2007). Values for 
NGC 1961 are shown in black for comparison. 


corrected circular velocity from the HI gas (Figure 10) by 
%/3. In this case, NGC 1961 falls squarely within the relation 
for early-type galaxies. 

We can show this in more detail with a detailed com¬ 
parison between the properties of NGC 1961 and compara¬ 
ble elliptical galaxies (Figure 14). For this comparison we 
consider the galaxies NGC 57, NGC 4455, and NGC 7796, 


studied by O’Sullivan & Ponman (20041 and O’Sullivan et 
[aL]p007| ). They examine IC 1531 as well, but this galaxy 
has an AGN with a radio jet, and they cannot constrain 
the hot halo properties very precisely. These galaxies were 
selected as massive ellipticals which are not in any known 
group or cluster and which have no neighbors within 2 B- 


magnitudes of their own magnitude (O’Sullivan & Ponman 
|2004[ ). |0’Sullivan et al.] ( |2007[ ) note that the well-known el¬ 
liptical NGC 720 falls near the border of their isolation cri¬ 
teria, but seems to have a higher density of faint galaxies 
around it than the three galaxies in their sample. As far as 
we can tell, NGC 1961 fulfills their isolation criteria as well, 
although it is less surprising for spiral galaxies to be isolated 
as compared to ellipticals. 

Within the uncertainties, the mass profiles and the en¬ 
tropy profiles for these four galaxies are all roughly consis¬ 
tent with one another. NGC 7796 has a slightly lower in¬ 
ferred total mass, so O’Sulivan et al. (2007) speculate that 
the hot gas around this galaxy might not be in hydrostatic 
equilibrium. There is a wider range of temperatures and 
abundances between these galaxies; NGC 1961 lies in the 
middle of the range of temperatures but has a lower abun¬ 
dance than the elliptical galaxies. 

An interpretation of the discrepancy with the Lx — Lk 
relation and the agreement with the Lx — cr relation is that 
the halo mass governs the X-ray properties of the system, 
instead of the stellar mass. As discussed in the Introduction, 
massive blue galaxies generally seem to lie in less massive 
halos than their red counterparts. As an example of this, 
from Anderson et al. (20151 we find that the effective halo 
mass of an average galaxy with the stellar mass of NGC 1961 


is around 6 x 10^^ Mq - much higher than the halo mass we 
measure for NGC 1961. The lower halo mass offers a natural 
explanation for why NGC 1961 appears to scatter so low on 
the Lx — Lk relation when compared to elliptical galaxies 
in much more massive halos. 


7 CONCLUSIONS 


In this paper we have analyzed the deepest X-ray observa¬ 
tion to date of a hot gaseous halo around a spiral galaxy. 
NGC 1961 has already been studied with both Chandra 
and XMM-Newton, but our deep observations allow us to 
perform a much more precise and robust analysis. Using a 


spatial analysis similar to Anderson & Bregman (20111 and 
Bogdan et al. ( 2013| ), we detect extended soft X-ray emis¬ 
sion to very large radii (at least 80 kpc). No breaks or bright 
features are evident in the surface brightness profile or the 
soft X-ray image of the galaxy. 

We also perform a spectral analysis on the XMM- 
Newton data. This allows us to measure the physical prop¬ 
erties of the hot gas out to 42 kpc. We explore eight differ¬ 
ent sets of spectral models, and the scatter between results 
from these different models allows us to estimate the sys¬ 
tematic uncertainty from the choice of models. In general 
this is comparable to the statistical uncertainties. The spec¬ 
tral and spatial methods give fully consistent results as well, 
further underscoring the robustness of our results. 

The hot halo has a negative temperature gradient (de¬ 
clining with radius), as is typical for elliptical galaxies of 
similar mass. The metallcity of the halo is consistent with 
being fiat at a value around O. 2 Z 0 . Such a low metallicity is 
not typical for elliptical galaxies of similar mass, which usu¬ 
ally have abundances around Solar, although X-ray fainter 
ellipticals can have sub-Solar abundances as well. The higher 
abundances for ellipticals are thought to stem from the mix¬ 
ing of the hot interstellar medium and the hot halo, whereas 
NGC 1961 can lock up the majority of its metals in the cooler 
phases of its ISM. In support of this picture, we find a sta¬ 
tistical preference for a two-temperature medium within the 
disk of this galaxy, with the dominant component having 
the high temperature and low abundance characteristic of 
the hot gaseous halo and a less massive component having 
a lower temperature and a higher abundance typical of the 
hot interstellar medium. 

We also measure deprojected density profiles for the hot 
halo out to 42 kpc. The results are consistent with the para¬ 
metric fit obtained by B13, when rescaled to our fiducial 
metallicity of O. 2 Z 0 , as well as with the much more uncer¬ 


tain profile obtained by Anderson & Bregman (20111. No 


evidence of a break or a flattening in the profile is seen. 
Extrapolating the profile to the virial radius, we estimate a 
hot halo mass which is comparable to the stellar mass of the 
galaxy. 

From the density profile and the temperature profile, we 
compute a pressure profile for the hot halo and then we in¬ 
fer a total mass profile under the assumption of hydrostatic 
equilibrium for the hot gas. This mass profile is roughly con¬ 
sistent with measurements from the circular velocity profile 
observed with 21 cm emission. Based on the extrapolated 
hot gas mass and total mass, the problem of missing baryons 
from NGG 1961 seems to persist, with an estimated baryon 
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fraction (hot gas + stars + neutral Hydrogen) around 30% 
of the Cosmic fraction. 

It is conceivable that photoionized gas at intermediate 
temperatures could also exist around this galaxy, as obser¬ 
vations by the COS instrument on HST suggest for many 
smaller galaxies. This could supply the missing baryons, but 
there are a few complications. First, such clouds would be 
far out of pressure equilibrium with the hot halo out to sev¬ 
eral hundred kpc. Second, the hot halo appears to be stable 
against global cooling instabilities over the range we can 
measure, making it difficult to condense such clouds out 
of the hot halo. These are not insurmountable difficulties, 
but they do complicate the modeling of the circumgalactic 
medium. 

Finally, we show that the X-ray luminosity from NGC 
1961 is about an order of magnitude fainter than one would 
expect from the Lx-Lk relation for elliptical galaxies. On 
the other hand, this galaxy seems to have exactly the right 
Lx for its dark matter mass, suggesting the dark matter halo 
governs the hot halo more fundamentally than the galaxy at 
the center. The non-detection of this galaxy in the Planck 
SZ map also supports this conclusion, as an elliptical galaxy 
with the same stellar mass and distance is very likely to have 
been detected. 

A number of open issues remain. It is critically impor¬ 
tant to push measurements out to larger radii, which will 
require deeper X-ray observations, possibly combined with 
an SZ and/or absorption-line analysis (the latter two are 
more sensitive at larger radii than X-ray emission, which 
depends on the square of the gas density). It is also impor¬ 
tant to increase the sample size of hot gaseous halos, in order 
to understand better how the properties of the gaseous halo 
connect with the properties of the galaxy and of the dark 
matter halo. 
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